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We study a three dimensional Z(3)-symmetric effective theory of high temperature QCD. The 
exact lattice-continuum relations, needed in order to perform lattice simulations with physical 
parameters, are computed to order 6[cP) in lattice perturbation theory. Lattice simulations are 
performed to determine the phase structure of a subset of the parameter space. 
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1. Introduction 

At high temperature, QCD matter undergoes a deconfinement transition, where ordinary had- 
ronic matter transforms into strongly interacting quark-gluon plasma. In the absence of quarks, 
Nf = 0, the transition is a symmetry-breaking first order transition, where the order parameter is 
the thermal Wilson line. The non-zero expectation value of the Wilson line signals the breaking of 
the Z(3) center symmetry of quarkless QCD at high temperatures. 

The transition has been studied extensively using lattice simulations [|I]], but becomes com- 
putationally exceedingly expensive at high temperatures T >5Tc. At high T, the complementary 
approach has been to construct perturbatively effective theories, such as EQCD and MQCD, using 
the method of dimensional reduction iQ]. In the dimensional reduction procedure, however, one 
expands the temporal gauge fields around one of the Z(3) vacua and thus explicitly violates the 
center symmetry and the effective theories fail to describe QCD for T below 57^. 

As a unification of these strategies, a 3D effective theory of hot QCD respecting the Z(3) 
symmetry has been constructed in [^]. At high temperatures, the effective theory is matched to 
EQCD but still preserves the center symmetry. The effective theory is further connected to full 
QCD by matching the domain wall profile separating two different Z(3) minima. Thus, one hopes 
that the range of validity of this theory would extend down to 7^. 

In order to perform lattice simulations the theory has been formulated on a lattice and the 
lattice theory has been matched to the continuum theory in [||]. The effective theory is super- 
renormalizable, and thus the exact relations between the continuum MS and lattice regulated theo- 
ries can be obtained via two-loop lattice perturbation theory. 

2. Theory 

The theory we are studying is defined by a three dimensional continuum action, which we 
renormalize in the MS scheme 

S = J d'-^'x I ^TrF^j + Tr {DiZ'^DiZ) + VoiZ) + Vi (Z) | , (2.1) 

where 

Fij = diAj-djAi + ig3[Ai,Aj], Di = di-ig3[Ai, ] (2.2) 

and Z is a 3 X 3 complex matrix, which in the limit £ ^ has dimension dimZ = \/GeV. The 
gauge fields A, are Hermitean traceless 3x3 matrices and can be expressed using generators of 
SU(3), Ai = A'fT", with Ti:T"T'' = ^5"^ . The covariant derivative is in the adjoint representation. 
The potentials are 

Vo{Z) = ciTr[Z+Z] +2c2Re(Det[Z]) +C3Tr[(Z^Z)2], (2.3) 

Vi{Z) = diTr[M'^M] + 2d2Re{Tr[M^]) + d3Tr[{M^M)^], (2.4) 

where M = Z — ^Tr[Z]l is the traceless part of Z. Here, the gauge coupling ^3 has a positive 
mass dimension dim[g3] =GeV, making the theory super-renormalizable. Because of the super- 
renormalizability, the coefficients 02,03,^2. and di, are renormalization scale independent and only 
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the mass terms ci and di acquire a scale dependence in the MS renormaUzation scheme. The 
coefficients c,, di, and §3 can be matched to the parameters of full thermal QCD by imposing 
the conditions that the theory reduce to EQCD at the high temperature limit, and that the theory 
reproduce the correct domain wall profile of full QCD [^]. This defines a subset of parameter values 
(with a limited accuracy due to perturbative matching), for which the theory describes thermal 
QCD. However, we consider here the theory in general, and do not restrict ourselves only to the 
perturbative matching regime. 

The action is invariant under local gauge transformations, with Z transforming in the adjoint 
representation: 



Aiix) 
Z{x) 



G{x)[Ai{x)--di]G-\x) 



■G{x)Z{x)G^\x 



(2.5) 
(2.6) 



where G{x) SSU(3). In addition to the local transformations, the action is invariant under the global 
Z(3) transformation 



Jlnn/^r^ 



1,2,.. 



(2.7) 



3. Lattice action 



The lattice action corresponding to the continuum theory can be written as S = Sw + 5z, where 



x,i<j 



1 - ^ReTr[^^v] 



(3.1) 



is the standard the Wilson action with the lattice coupling constant j3 = 6/{ag'^). 

The kinetic term, Tr {DiZ^DiZ) , is discretized by replacing the covariant derivatives by covari- 
ant lattice differences. Then the scalar sector of the action reads: 



5z = 2 ( -^ I £ReTr \z''Z-Z\x)Ui{x)Z{x + T)U^ {x) 



(3.2) 



+ ( "ft ) L (ciTr[Z'^Z] + 2c2ReDetZ + c^Tv[{Z^Zf] + diTv[M^M] + 2J2ReTrM^ + d^Tv[{M^Mf]) 

where Ci,di,M, and Z are dimensionless. Only the mass terms ci and di require non-trivial renor- 
malization and all the other terms can be matched to order 0{(P) on tree-level by simply scaling 
withes: 



Z = g^Z, M = g^M 

C2=glc2, d2=gld2 
C3 =§3^3, d3 =gld3. 



(3.3) 
(3.4) 
(3.5) 



For the mass terms, renormalization has to be carried out, so that the long distance physics is the 
same in both regularization schemes. A two-loop lattice perturbation theory calculation gives: 



^■=i-s"5'«^^«^''^-i^ 



64c3 + —ci ] (log j8 +0.08849) + 37.0863C3 



+ 



-In 



(3.6) 
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and 



=4- — (3-17591 + 5.64606^4) ^(41.780852 + 37.0863^4 



( —cj - 64^3 + -^^3^3 + y ^1 + l) [logj8 + 0.08849] ] + &{p-'). (3.7) 



There are also higher order corrections (corrections of order ^(j8^^) corresponding to order ff{a) 
in lattice spacing), but their effect vanishes in the continuum limit. Various operators have also 
been renormalized in [Q] on the lattice in order to convert their expectation values to continuum 
regularization. 

4. Phase diagram of Vi (Z) 

A simpler model is obtained from the original theory by setting c; = 0. In this model, the trace 
of Z decouples and can be integrated over as a free scalar field and the relevant degree of freedom 
is thus a traceless complex matrix M. The model is defined by the action: 



d^x 



-TrFJ^j + TrDiM^DiM + diTrM^M + 2J2Re(Tr[M^] ) + d3Tr{M'^Mf 



(4.1) 



If the cubic term d2 is zero, the Lagrangian is invariant under a U(l) global symmetry M —>■ gM, 
g SU(1). The breaking of the symmetry is signalled by a local order parameter: 



(Tr(M + Mt)3)2 + (Tr(M-Mt)3)2. (4.2) 

In the symmetric phase £/ is strictly zero and in the broken phase it has a non-zero vacuum ex- 
pectation value, while the two phases are separated by a first order transition. In the broken phase 
(TrM^M) is larger than in the symmetric phase. After the inclusion of the cubic term, £/ is no 
longer strictly an order parameter, since the U(l) symmetry is explicitly broken. However, the first 
order transition remains and is accompanied with a significant discontinuity in £/ and (TrM ^M) . 

4.1 Lattice analysis 

A non-perturbative lattice analysis has been performed to obtain the phase structure of the 
model. For the simulations we used a hybrid Monte-Carlo algorithm for the scalar fields and 
Kennedy-Pendleton quasi heat bath and full group overrelaxation for the link variables. 

The transition was found to be of the first order for all parameter values used in the simulations 
(ds < 4 and fi?2 < 0.15) accompanied with a large latent heat and surface tension; hysteresis curves 
showing discontinuity around critical point in (TrM^M)-^ can be seen in Fig.|I| The probability 
distributions of TrM^M along the critical curve are very strongly separated (see Fig.^. This makes 
the system change its phase very infrequently during a simulation, and a multicanonical algorithm 
is needed to accommodate a phase flip in reasonable times for any system of a modest size. Even 
with the multicanonical algorithm, the critical slowing restricts us to physical volumes up to V < 
50/gf. 

The pseudo-critical point was determined requiring equal probability weight for TiM^M in 
both phases. The simulations were performed with j8 = 12 and a lattice size N^ = 12^, which 
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Figure 1: Discontinuity in the quadratic condensate in continuum regularization {TrM^' M)-^ for d^ = 
0.1,1,3. The phase transition gets weaker as the coupling ^3 grows. The metastable regions shrink and 
the discontinuity diminishes. 



precludes the continuum extrapolation as well as the thermodynamical limit. However, these limits 
were studied for one set of parameter values. The dependence of the critical point on the lattice 
spacing was beyond our resolution for the lattice spacings used and the volume dependence was 
found to be of order of five per cent for the volumes used (see Fig.|| ). 

The phase diagram can be seen in Fig.^ The non-perturbative critical line measured from the 
lattice follows the one-loop perturbative result for small values of ^3, but for larger di, fluctuations 
make the system prefer the symmetric phase. The discontinuity in (TrM^M) along the critical 
line diminishes, as dj, gets larger, but it seems that the discontinuity persists, even if its magnitude 
diminishes in the Umit d^ ^0° suggesting that there is a first order phase transition for any (positive) 
value of J3. 




Tr(MV)/g^ 



Figure 2: Histograms of TrM^M in logarithmic scale with d2=0 along the critical curve. Transition channel 
between the peaks weakens and the transition gets stronger for decreasing 1^3. For dj,/g^ = 0.5, the relative 
probability density in the tunneling channel is suppressed by a factor ^ 10^'". 



5. Conclusions 



The exact relations between the lattice and continuum MS regulated formulations of the Z(3)- 
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Figure 3: Volume and lattice spacing dependence of the pseudo-critical point with dj,=2 and ii2 = 0. 1 . The 
pseudo-critical point was determined by requiring equal probability weight for TiM^M in both phases. The 
lines represent linear fits. The dependence on lattice spacing and volume seem to be within 5% for the lattice 
spacings and volumes used. 



symmetric 3D effective theory of hot QCD have been calculated in |Q]. The Lagrangians and the 
operators up to cubic ones have been matched to ^(a"). These results make the non-perturbative 
lattice study of the theory possible. 

An interesting model with non-trivial dynamics is obtained by setting c, = in Eq.(2.3). The 




Figure 4: The phase diagram of the soft potential as a function of di,d2 and d^. First order critical line 
separates two phases. Solid lines represent polynomial fits to the lattice data points and dashed lines are the 
perturbative predictions. The symmetric phase refers to the phase where with d2 — the order parameter 
vanishes and with d2^0is smaller than in the broken phase. On the right panel d$ — 2. 
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phase diagram of this model has been determined using lattice simulations. Two distinct phases 
were found, separated by a strong first order transition. 

In the future, it is our goal to map out the phase diagram in the full parameter space of the 
theory, in order to search for regions in which the phase diagram would resemble that expected for 
the finite-temperature SU(3) pure Yang-Mills theory. 
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